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Abstract: In this paper we complete our re-assessment of the production of W boson 
pairs at the LHC, by calculating analytic results for the gg — > W + W~ — > vl + l~v process 
including the effect of massive quarks circulating in the loop. Together with the one-loop 
amplitudes containing the first two generations of massless quarks propagating in the loop, 

q^ ■ these diagrams can give a significant contribution with a large flux of gluons. One of the 

component parts of this calculation is the production of a standard model Higgs boson, 
gg —7- H and its subsequent decay, H —> W + (— > vl + )W~ {— > l~v). We will quantify 

i^ the importance of the interference between the Higgs boson production process and the 

gluon-induced continuum production in the context of searches for the Higgs boson at the 
Tevatron and the LHC. For instance, for run < 140 GeV the effect of the interference 
typically results in around a 10% reduction in the expected number of Higgs signal events. 
The majority of this interference is due to non-resonant contributions. Therefore cuts on 
the transverse mass such as those currently used by the ATLAS collaboration reduce the 
destructive interference to about a 1% effect. We advocate that a cut on the maximum 
transverse mass be used in future Higgs searches in this channel. 
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1. Introduction 

The search for the Standard Model (SM) Higgs boson is entering the closing stages. With 
the LHC acquiring data at an impressive rate, evidence for the existence of the Higgs boson 
or its exclusion over the mass range 115 — 600 GeV should be expected by the end of 2012. 
In addition, ongoing studies at the Tevatron with around 10 fb~ will provide additional 
independent exclusion regions over a wide range of potential Higgs masses. In summary, 
the combined results from both hadron colliders can be expected to tightly constrain the 
Higgs boson within the course of the coming year. Winter 2010 results from the various 
collaborations can be found in refs. [1—3]. 

An important Higgs search channel over the range 120 — 200 GeV is the process 
gg — >• H — >• WW — > vU'v' . This proceeds through gluon fusion via a top (or bottom) 
loop with the subsequent production of two W bosons that then decay leptonically. In 



setting a limit on this channel one needs to possess precise predictions for the cross section 
and kinematic distributions. Thankfully, as a result of much theoretical work the (fully 
differential) production cross section for a Higgs boson through gluon fusion is known at 
NNLO (see for example refs. [4-10]). By including the NNLO corrections one typically 
finds scale uncertainties on inclusive quantities (which are used to quantify theoretical un- 
certainty by the collaborations) of O(10 — 20%) [11]. Quoting an error of this order requires 
that all contributions to the cross section that could possibly change it by 0(10%) have 
been calculated, and included in the theoretical predictions used. To correctly exclude 
the Higgs then, one must confirm that no remaining O(10%) effects are neglected in the 
theoretical prediction. 

One potentially large contribution to the cross section is the interference between the 
Higgs signal and SM continuum production of WW pairs. Power counting in as indicates 
that this contribution is the same order as the LO Higgs signal cross section. Typically these 
interference terms are neglected in the theoretical calculations used by the collaborations. 
The aim of this paper is to fully quantify this interference under a variety of experimental 
cuts. In order to do this we calculate analytic results for the gg —> WW — > vl + l~D process 
including the effect of massive quarks circulating in the loop. The interference has been 
studied in ref. [12], in which the focus was on a 14 TeV LHC, and on the effects on the 
total cross section. In this study we will re-address the issue in the context of current Higgs 
boson searches. The interference effects associated with Higgs production and decay to ZZ 
have been studied some time ago [13, 14]. In these papers the primary focus was at SSC 
operating energies and large Higgs masses. Similar studies have also been performed for 
the case of a light Higgs boson decaying into di-photons at the LHC [15] and the case of 
Higgs production at a photon collider with subsequent decay into bottom quark pairs [16]. 

The process qq — > WW was first calculated in the Born approximation in ref. [17], 
with strong corrections later computed in refs. [18-21]. Phenomenological NLO results 
for the Tevatron and the LHC have been presented in refs. [22-24]. The gluon-initiated 
contribution gg — »• WW was first calculated in refs. [13, 25] and later improved to account 
for off-shell effects of the vector bosons and their subsequent decays in ref. [26] . The effect 
of massive quarks circulating in the loop was later assessed in ref. [12]. Combined results, 
including both the qq and gg initiated contributions - for loops of massless quarks only - 
were recently presented in ref. [24]. In this paper we shall use modern methods to provide 
a calculation of the massive quark contributions that is equivalent to the calculation in 
ref. [12] and update the results of ref. [24] accordingly. 

The plan of this paper is as follows. In section 2 we shall describe the calculation 
of the amplitude gg — > WW where a massive top quark (and, where numerically signif- 
icant, a massive bottom quark) circulates in the loop. Sections 3, 4 and 5 will present 
phenomenological results using the implementation of these amplitudes into the parton 
level Monte Carlo program MCFM. Section 3 will address the size of the corrections to the 
WW process resulting from the inclusion of massive quarks. In sections 4 and 5 we shall 
discuss the interference contributions at the LHC and Tevatron respectively. In particular 
we will discuss these effects in the region where cuts are applied to enhance the signal to 
background ratio for the detection of the Standard Model Higgs boson. In section 6 we 



discuss the impact of the interference in regards to the NNLO predictions for the total 
cross section. Finally, in section 7 we draw our conclusions. 

2. Calculation 

In reference [24] we showed that the calculation of gg —> WW with massless quarks circu- 
lating in the loop can be extracted from Ref. [27]. Specifically, we find that the contribution 
from a single generation of massless quarks in the loop is given by, 
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The helicities and colour labels of the two gluons are hi, hi and a±, a-i respectively and the 
propagator factor is, 

V w (s) = ^— , (2.2) 

where m\y and T\y are the mass and width of the W boson. The amplitude Ag. 4 and 
constant cr are defined in Sections 2, 6 and 11 of ref. [27]. The particle labelling on the 
right hand side of this equation is as written in ref. [27]. For our purposes in Eq. ( |2.l| ) we 
make the identification, (q — >• £, q — >• v#, e — )■ £' and e — > Vgi) in going from right to left. 

The calculation in Ref. [27] was performed using unitarity methods. Our aim in this 
section is to perform a similar calculation but including the mass of the top quark circu- 
lating in the loop and also the mass of the bottom quark where numerically significant. 

2.1 Amplitudes for gg — > WW 

In this section we present results for the amplitudes relevant for the process, 

-> g(pi) + g(p2) + MPs) + £{pi) + i'ips) + %(pe) • (2.3) 

Topologies of diagrams that could potentially contribute to the leptonic final state in which 
we are interested are shown in Fig. ffl. When the vector boson couples to the triangular 
fermion loop (diagrams (a) and (b)) via a vector coupling, the contribution vanishes via 
Furry's theorem. When the vector boson couples via an axial coupling to massless quarks, 
the contribution of up-type quarks with weak isospin 1/2 and down-type quarks with weak 
isospin -1/2 cancel. In the case of the third generation with a massive top quark this 
cancellation is no longer operative, but it turns out that there is no net contribution [12] 
as we shall now illustrate. If the masses of the appropriate lepton pairs are constrained to 
have the mass of the W, the contribution of diagram (a) vanishes and diagram (b) is not 
present. This follows because the only non- vanishing contribution of the triangle diagram 
is proportional to p% where pz is the momentum of the Z boson and \i is its Lorentz 
index [28]. If this condition is relaxed to allow off-shell VF-bosons, the contribution of an 
individual fermion flavour cancels between the doubly resonant diagrams (a) and the singly 
resonant diagrams (b). Therefore in both cases the triangle diagrams vanish because of 





Figure 1: Topologies of diagrams that could potentially contribute the process gg — > v e e + > v^. 



electroweak gauge invariance. Thus the only triangle loop contribution comes from the 
Higgs boson mediated diagram (c). 

We shall therefore separate the calculation of the full amplitude as follows, 
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where the first two terms in the square brackets represent the six continuum diagrams 
shown in Fig. ^. We explicitly separate out the contribution from the first two generations, 
(2^4 mass i ess ), in which we consider the quarks in the loop to be massless, and the contribu- 
tion from the third generation in which the top quark mass is not neglected, -Amassive- We 
have also extracted an overall color, coupling and propagator factor that is the same as in 
the massless case, so that comparing with Eq. ( |2.lD we immediately see that, 
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The final contribution .Amggs originates from two triangle diagrams of the type shown in 
Fig. |(c). 

2.2 Integral basis 

We consider the one-loop amplitude -Amassive, corresponding to third generation quarks 
circulating in the loop. In our calculation we retain the dependence on the mass of the top 
quark but treat the 6-quark as massless since the effect of a non-zero value for m& is at the 
level of 0.1% [12]. 

The amplitude can be expanded in terms of scalar integrals and a rational part as 
follows, 
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Figure 2: Diagrams that contribute to gg — > WW with the third generation of quarks running in 
the loop. The top quark is denoted by a heavy (magenta) line whereas the b bottom quark is shown 
in blue. The same diagram topologies contribute to the case where first and second generation 
quarks run in the loop. 
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Table 1: Definitions of the scalar integrals that appear in the calculation of the amplitude A„ 

i.e. continuum production of gg — > W + W~ through a loop containing at least one massive particle. 



The scalar integrals BW, C®, D^> that contribute to this particular amplitude are listed in 
Table 111. The coefficients &,, c,, di are amenable to calculation by unitarity cuts; the contri- 
bution 1Z represents the remaining non-cut-constructible rational part. We have performed 
a complete independent numerical check of our analytic calculation of the coefficients and 
the rational part using D-dimensional unitarity [29]. No tadpole contributions are found 
to be present. 



The D-dimensional scalar integrals themselves are defined as follows, 
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where g n = X^iLi^ anc ^ ^^i = (Pi+Pj) 2 - ^ or t^ 16 purposes of this paper we take the masses 
in the propagators to be real. Near four dimensions we use D = 4 — 2e (and for clarity 
the small imaginary part which fixes the analytic continuations is specified by +ie). /x is a 
scale introduced so that the integrals preserve their natural dimensions, despite excursions 
away from D = 4. We have removed the overall constant which occurs in D-dimensional 
integrals 
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The final numerical evaluation of the amplitudes uses the QCDLoop library [30] to provide 
values for these scalar integrals. 

2.3 Rational part 

In calculating the box diagrams we shall start with the least challenging terms in Fig. 



using the decomposition of Eq. (2.6). The rational terms are completely determined by the 
massless calculation which has been already presented. The box diagrams are of rank four, 
i.e. they contain at most four powers of the loop momentum in the numerator. However 
the first mass dependent terms contain at least two powers of the mass from numerator 
factors and are therefore at most of rank two. These contributions are cut-constructible 
[27] and do not contribute to the rational part. The rational terms are controlled by the 
ultraviolet behavior and thus the presence of masses in the denominators of the rank 3 and 
4 boxes does not change the rational part. The rational part is therefore identical to the 
massless case. This assertion has been checked by direct numerical evaluation. 

2.4 Bubble coefficients 

Turning now to the coefficients of the bubble integrals, we find numerically that the results 
for the massless case can be re-used for the massive case. In the limit that the internal 
mass m vanishes, we see that the integrals B^ l > and B^ & > become degenerate. Therefore 
the integral B^ 6 > can be dropped from the basis and we can write the bubble contributions 



to the amplitude for the massless calculation as a sum over only 5 terms, 

massless massless ' v / 

i=l,5 

In this equation the subscript "massless" labels the coefficient in the massless calculation 
and ^massless = ^ \m=o- We find that four of the integral coefficients are identical to the 
massless case, 

^ = Csiess for i = 2, 3, 4, 5. (2.10) 

Moreover, the sum of the two remaining coefficients is equal to the final massless coefficient, 

as a result of the fact that the amplitude does not contain any poles of ultraviolet origin, 

Y^ b {i) = . (2.12) 

i=l,6 

We have performed an explicit calculation of the simplest coefficient (&' 1 ') (using the tech- 
nique described in [31]) and have hence determined all six of the bubble coefficients from 
the known massless bubble coefficients in Ref. [27]. 

2.5 Box coefficients 

We express the amplitudes in terms of spinor products defined as, 

(ij) = u_{pi)u + {pj), [ij] = u + {pi)u_(pj), (ij) [ji] = 2p, L ■ pj , (2.13) 

and we further define the spinor sandwich, 

(i\(j + k)\l] = (ij)[jl] + (ik)[kl}. (2.14) 

In order to express the amplitudes compactly we define two symmetry operations, 

fli Pl := {(3*»6),(4*»5),<>*»Q}, (2.15) 

flip 2 := {(1 O 2), (3 ++ 6), (4 o 5), () O []} , (2.16) 

which we will use to relate various coefficients in the amplitude. Specifically for the box 
coefficients in Eq. ( |2.6[) related to the box integrals with the labels given in Table [I] we 
have, 

d ti = fli Pi d i,5 ' d 2,e~ = fli Pi d ti' d 2,i = ffipi d tfT' d i,6 = fljpi d i,i ■ ( 2 - i7 ) 

Finally the relations between the coefficients for boxes 3 and 4 are, 

4+ = flipa d 3 -, d 4 - = flip 2 4+, dj + = flip 2 dz + , d\~ = flip 2 4" . (2.18) 

Note, in particular, that the flip 2 symmetry relates the (— ,+) coefficients of boxes 3 and 
4 whereas flip 1 relates d 2 + 1° d \ • 



Given these relations, it is therefore sufficient to compute the coefficients d\,d^ and 
d$, which we do using the quadruple cut method [32]. The simplest coefficients are those 
of the boxes containing only a single massive propagator, 

4 + = , 
d~ = , 

<k ="2^^ (1|(3 + 4)|2]4 ((13) [34] (56) [26] +m (15) [42]) 
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2 S34S56 (2|(3 + 4)|lj 4 
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The coefficients of the boxes containing two massive propagators are more complicated but 
still relatively compact, 
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The coefficients of box 5 are lengthier; interested readers may inspect their form in the 
distributed MCFM code. 

2.6 Triangle coefficients 

Turning now to the coefficients of triangle integrals, we see from Table [l] that twelve 
coefficients are required. These triangle integrals are depicted in Figure |3[ Since the 
amplitude that we are calculating is finite, all of the poles in e must cancel. This cancellation 
of infrared poles leads to relations between the triangle coefficients and the singular boxes 
1-4. In particular we find that c±, c 2 , C3, C4 and C5 can be determined in terms of these 
box coefficients through the set of equations, 

-^- 2 +^^- 2 + Cl = 0, (2.21) 
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Figure 3: Scalar triangle integrals present in the calculation. 
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(2.25) 



A ' S12(S156 -m 2 ) S156-S34 

These relations involve a kinematic factor A related to boxes 3 and 4, where A is given by, 

A = (m 2 - s 13i )(m 2 - si 56 ) - (m 2 - s 34 )(m 2 - s 56 ) . (2.26) 

Therefore it only remains to determine c« for i = 6-12. Many of these coefficients are related 



to one another by the flip operations defined in Eqs. ( 2.15 , 2.16|) or simpler symmetries. In 
particular we find that the different helicity amplitudes for triangle 6 are related, 
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c+~(lf*2) 



(2.27) 



whereas the flipj^ symmetry of Eq. (2.15) relates the coefficients of triangles 7 and 9, 
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fliPi C++ , 
fliPi Cj~ , 
flipx c+~ , 



c+~ = fli Pl cj + . (2.28) 

The coefficients cs and c\q can be determined from those for cj and eg using the relations, 
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(2.29) 

Finally, we have the relationships between the different helicity amplitudes for triangle 11, 

c u - = c++{(l o 2), (3 *» 4), (5 *» 6),() «■ []} , 

c n + = c+f(lo2), (2.30) 

and triangle 12, 

C\2 = fli Pi C 12 + > 

c u + = fh Pl c+- . (2.31) 

We have calculated the coefficients of the remaining unconstrained triangles using the 
method of [33]. The results for these triangle coefficients are complicated but may be 
inspected in the distributed MCFM code. 

2.7 Calculation of the Higgs mediated diagrams 

The calculation of Higgs boson production via a loop of heavy quarks is well known. Here 
we also include the decay of the Higgs boson into a pair of W bosons, with our usual 
particle labelling. The contribution to the full amplitude is then, 

<i g gs = K hlM ViiM E [ C #K) h(s 12 ,0,0)m 2 q ,m 2 q ,m 2 q )+n H (m q j\ , (2.32) 

q=t,b 



where overall factors of color, couplings and W propagators have been extracted in Eq. (2.4). 
Note that here we sum over both top and bottom quark loops since the effect of including 
the bottom quark contribution is numerically important for lower Higgs masses (through 
its interference with the top quark loops). Although we have therefore treated the bot- 
tom quark mass in an inconsistent manner (finite in the loops coupling to a Higgs boson, 
zero otherwise), we have checked that the discrepancy is not phenomenologically relevant. 
Although including the mass of the bottom quark can significantly change the total pro- 
duction rate for a light Higgs, the interference terms that are the focus of this paper have 
a very small dependence on nib. The dependence on the helicities of the gluons is encoded 
in the factor K u 2 which is given by, 

K~+ = K+- = , 

2(12) (3 5) [64] 



K 
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The coefficient of the scalar triangle and the rational contribution are, 

. . im 2 ( Am 2 \ „ . . im 2 
cn{m) = — 1 , TZ H {m) = . (2.34) 
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Substituting for /3(si2,0, 0;m 2 ,?7i 2 ,m 2 ) in Eq. ( 2.32 ) we obtain, 
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where the function F(x) is defined by, 
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-2[sin- 1 (l/2Vi)] : 



x < 1/4, 
x > 1/4. 



(2.36) 



Up to overall factors that are associated with the decay H — > WW, this is the well-known 
result (see e.g. Ref. [34]). 

3. Effect of massive loops in gg — > WW 

3.1 Parameters 

The results presented in this paper are obtained with the latest version of the MCFM code 
(v6.1). We use the default set of electroweak parameters that assumes the following set of 
inputs, 



m w = 80.398 GeV , m z = 91.1876 GeV 
IV = 2.1054 GeV , T z = 2.4952 GeV , 
G F = 1.16639 x 10~ 5 GeV~ 2 . 



(3.1) 



Using the values of m^y, mz and Gf as above then determines a e , m ,(mz) and sin 2 9 W as 
outputs, where 9 W is the Weinberg angle. We find, 



sin 



a,. 



.(rn z ) 



1 - m^/mz = 0.222646 , 
V^GFrn^y sin 2 8 W 1 



IT 



(3.2) 



132.338 

Where massive loops of top and bottom quarks are included we use nit = 172.5 GeV and 
m.6 = 4.4 GeV. When we include the contribution from diagrams involving a Higgs boson, 
the values of the width that we use are taken from HDECAY [35] and are shown in Table |2|. 
For the parton distribution functions (pdfs) we use the sets of Martin, Stirling, Thorne 
and Watt [36]. We use the NLO pdf fit, with a s (mz) = 0.12018 and 2-loop running. In 
this section we will use a common renormalization and factorization scale equal to mw- 
In the following sections, where the interference with the Higgs amplitudes is studied, the 
common scale (unless otherwise stated) will be set equal to mu- 

3.2 Results 

We begin by assessing the numerical impact of the contribution arising from the loops con- 
taining massive quarks, i.e. the (£, b) doublet represented by the term ^4 maS sive m Eq. (J2~i|). 
This contribution was neglected in the analysis of ref. [24] based on earlier studies at the 
14 TeV LHC where it was found to be very small [12]. 
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m H [GeV] 


T H [GeV] 


m H [GeV] 


T H [GeV] 


100 


0.2573 xl0~ 2 


220 


2.301 


110 


0.2938xl0~ 2 


240 


3.397 


120 


0.3600xl0- 2 


260 


4.767 


130 


0.5006xl0~ 2 


280 


6.443 


140 


0.8281xl0~ 2 


300 


8.452 


150 


0.1744X10- 1 


320 


10.81 


155 


0.3042X10" 1 


340 


13.50 


160 


0.8255X10" 1 


360 


17.57 


165 


0.2434 


380 


23.04 


170 


0.3760 


400 


29.16 


180 


0.6291 


450 


46.82 


190 


1.036 


500 


67.94 


200 


1.426 


600 


122.5 



Table 2: Total Higgs width as a function of the Higgs mass using the HDECAY code (v3.51) from 
ref. [351. 



For now we neglect the contribution .Amggs in Eq. ( |2.4| ) and consider only the effect 
arising from adding the amplitude -4 mass i V e to the result for two massless doublets of quarks 
only. We note that adding the massive doublet results in new contributions from both 
"4-massive and the interference term, 2Re (^massive-^massiess)- ^ n the resm ts that we present 
here we will compare three contributions to the WW final state: 



a gg [n gcn = 2]: the gluon- initiated contribution for two massless doublets 



resulting from the term 2„4. mass i ess in Eq. (2.4); 



" gg l^gen 



3]: the gluon- initiated contribution for all three doublets 



resulting from the combination 2„4 mass i ess + A v 



in Eq. (I 



a tot"° : ^e total prediction obtained by adding a gg [n gen = 3] and 
the NLO result (qq, gq, gq initial states) [22, 24]. 

Our results for these cross sections, for the production of the final state W + {— > 
v e e + )W~ {— > fi^u^), are shown in Table |3|. We consider various values of the hadronic 
collider energy y/s, with proton-antiproton collisions at y/s = 1.96 TeV and proton-proton 
collisions otherwise. In these calculations we have applied no cuts to the final state. 

We first observe that including the third generation has only a minor effect on the 
gluon-gluon initiated contribution, ranging from an increase of 6% at the Tevatron to an 
increase of 13% at the 14 TeV LHC. This small increase is washed out even further when 
considering the correction to the total cross section including qq,gq and gq contributions. 
As expected, the gluon-gluon contribution is negligible at the Tevatron (0.4% of the total) 
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V^ [TeV] 


1.96 (pp) 


7 


8 


10 


12 


14 


@gg l^gen — ^J 
" r gg l^gen "J 
® gg L^gen "J/ "gg L^gen ^J 


0.460(0) 

0.490(1) 

1.065 


13.74(1) 

15.16(1) 

1.103 


18.19(1) 

20.12(2) 

1.106 


28.37(2) 

31.61(3) 

1.114 


40.06(3) 

44.84(4) 

1.119 


52.99(4) 

59.59(4) 

1.125 


CggPgen = 3]/<7 tot 


134.6(2) 
0.0036 


539(1) 
0.028 


657(1) 
0.030 


904(1) 
0.035 


1162(1) 
0.039 


1429(2) 
0.042 



Table 3: Cross sections for the process W + (-^ v e e + )W~(-^ ix~v^) (in femtobarns) with no cuts 
applied. The cross section resulting from the gluon-gluon initial state (a gg ) is broken down by the 
number of generations of quarks circulating in the loop. The result with three generations is also 
compared with the total cross section (a^°), obtained by also including the NLO corrections to 
the quark-antiquark initial state. 
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Figure 4: The ratio of a gg [n gen = 3] to a gg [n gen — 2] as a function of the top quark mass m t . The 
expected results for two and three massless generations are shown as dashed blue (lower) and red 
(upper) lines respectively. 



and even at the LHC it is at most 4% of the total at the highest energy foreseen (although 
this can increase upon application of various cuts on the final state [12, 24]). We note that 
our findings at 14 TeV are in complete agreement with the results presented in Ref. [12]. 

The fact that the three-generation cross section is so close to the two-generation result 
may be somewhat unexpected. For that reason, in Fig. H we illustrate the dependence of 
Pgg[ngen = 3] on the top quark mass, m t . In the limit of a very light top quark the third 
generation is equivalent to the two massless generations, and as a result cr gg [n gen = 3] — > 
(Jgg[n gen = 2] x 9/4 as m t — > 0. In the opposite limit, when the top quark is very heavy, 
the third generation decouples so that a gg [n gen = 3] — > a gg [n gen = 2] as mt — > oo. The 
behaviour of the three-generation result changes rapidly in the region mt ~ mw, as could 
be expected from the kinematics of the final state. 

In Fig. [| we present the px distribution of one of the W bosons for the first two gener- 
ations (red) and third generation only (blue) . We observe that at low pt the contribution 
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Figure 5: The transverse momentum distribution of the W bosons at y/s — 7 TeV and with no 
cuts applied. The contribution from gg — > WW through a loop of third generation quarks only is 
shown in solid red and the contribution from gg —>■ WW through two massless generations is shown 
in dashed blue. The distribution obtained from the NLO calculation of the qq, gq, gq processes is 
shown as the dotted magenta curve. 



from the third generation is negligible but in the region px > rat it becomes increasingly 
more important. This indicates that for searches that require a high-p^ W boson, the 
top loops should not be neglected. Numerical instabilities exist in the massive (and to 
a lesser extent in the massless) calculation in the region of low p^ . The origin of these 
instabilities can be traced back to terms of the form 1/(2|3 + 4|l] n in the n-point integral 
coefficients, for instance in Eq. ( |2,19| ) (for n = 4). Such denominators clearly vanish in the 
limit that the momentum (3 + 4) is collinear to either particle 1 or 2, i.e. the momentum 
Pw = P3 +Pa is along the beam, so that pt(W) = 0. The amplitude may be re-expressed in 
terms of combinations of the scalar integrals that are finite in this limit, a refinement that 
has already been partially performed for the massless amplitudes [27]. We find that some 
numerical instability still remains for the massless amplitudes and, in the massive case, 
the denominators are both more prevalent and more cumbersome to remove. Therefore 
to ensure numerical stability we instead simply impose a cut, pjf > 2 (0.05) GeV in the 
massive (massless) calculation. The application of this cut changes the total cross section 
by at most 0.05%. 



4. Higgs —7- WW interference effects at the LHC 

Since we have seen that that the effect of the gluon-gluon initial state is largest at the 
LHC we shall first present phenomenological results applicable for Higgs boson searches 
at yfs = 7 TeV. In particular in this section we will consider interference effects between 
terms in the amplitude that result from diagrams involving a Higgs boson and those that 
do not. 
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cut selection 


total cross section 


search cuts 


M H [GeV] 


140 


170 


200 


140 


170 


200 


on 


79.90 


116.2 


75.39 


1.884 


12.95 


1.664 


°H,i 


72.49 


114.4 


74.52 


1.758 


13.84 


1.999 


aH,i (extracted from Ref. [12]) 


72.56 


114.6 


74.45 


1.760 


13.87 


1.999 



Table 4: A comparison of results obtained from our calculation in MCFM with the corresponding 
values extracted from Table 3 of Binoth et al. [12]. Both W's decay leptonically into a single lepton 
flavor, no cuts are applied, y/s = 14 TeV and cross sections are given in fcmtobarns. The cross 
sections are computed either excluding (<jh) or including {(JH,i) the effect of interference with the 
gluon-initiated background process. 



The cross sections that we shall compute are obtained from Eq. (2A) by selecting 
various terms in the square of the amplitude, integrating over the final-state phase space 
and convoluting with parton density functions appropriately. We shall divide the full 
contribution to the cross section as follows, 



0\B 



\A 



box 



»^bo 



1A 



massless 



+ A 



massive i 



&H ► |-4ffiggs| , 

Gi — ► 2Re (^Higgs^box) 

&H,i = &H + <?i ■ 



(4.1) 



For a given Higgs mass the total cross section including the Higgs diagram is thus the 
sum of the first three terms, ob and on are the usual leading order gg background and 
Higgs signal cross sections respectively. When the interference is included, we consider the 
quantity OH,i as representing the cross section for events containing a Higgs boson, to be 
compared with the background-only cross section as- We note that, since an,i is not a 
physical cross section and therefore not constrained to be positive, we may (and indeed 
do) find that ajj,i < &h- 

To validate our computation of the interference terms we first compare our results with 
those presented in Ref. [12] for the LHC energy of y/s = 14 TeV. As already noted, in this 
reference the dependence on the bottom quark mass is retained but here we have chosen 
to set it to zero. The expected difference, of the order of 0.1%, is at the same level as the 
Monte Carlo uncertainty in our result. Within these uncertainties, the two calculations are 
in complete agreement for a range of Higgs masses and cuts, as can be seen from Table |3|. 
For the purposes of this comparison all parameters in MCFM have been set to the values 
used in Ref. [12]. Ref. [12] also specifies the Higgs search cuts used in the final three 
columns of the table. Although the quantity au,i is not explicitly given in Ref. [12] it may 
easily be extracted from the results presented there. 

4.1 No final state cuts 

We now return to the case at hand, namely the LHC operating at 7 TeV. The quantities 
(Xh and (JH,i, with no cuts applied on the final state, are shown as function of mjj in Fig. |(| 
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Figure 6: Upper panel: The cross sections for gg — > H — >• W + {-^ v e e + )W~(— > [i~v^) in femto- 
barns, with (oha) and without (cth) the interference with SM gg — » WW production. The dashed 
line represents the calculation of an,i including only the hrst two generations of quarks. Lower 
panel: The ratio of the cross sections with and without the interference terms. The dotted magenta 
line highlights the boundary between constructive and destructive interference. 



M H [GeV] 


120 


140 


170 


200 


400 


OH 


7.90(1) 


20.29(1) 


26.13(2) 


14.69(1) 


4.23(1) 


&H,i 


6.73(1) 


19.04(1) 


26.25(2) 


14.96(1) 


4.16(1) 




0.852 


0.938 


1.005 


1.018 


0.983 



Table 5: Cross sections for gg — > H — > W + (— > v e e + )W~ {— > [i^v^ in fcmtobarns at yfs = 7 TeV 
with no cuts applied, computed at leading order and either excluding (an) or including {<jh,i) the 
effect of interference with the gluon-initiated background process. 



Numerical values of these cross section are shown in Table |5| for a selection of benchmark 
Higgs masses. We observe that the relative size of the interference is strongly dependent 
on the Higgs mass and that the interference changes sign at the ran = 2mw threshold. 
For ran > 2mw there are two further changes of sign, with a minimum at m# = 2mt- For 
very large Higgs masses the interference becomes large and positive. For reference we have 
also plotted in Fig. |6| the contribution to the interference from the first two generations of 
quarks only (i.e. setting *4 ma ssive = in the definitions of Eq. (fP|). The difference between 
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the two and three generation results is striking. Except for values of mjj < 200 GeV, in 
general either the sign or the magnitude of the two-generation interference is changed in the 
three- generation result. This is a result of features in the three-generation result around 
mu = 2?7Jt, as is expected from the thresholds of the integral functions. We conclude that, 
despite being a negligible factor in total rates, it is imperative to include the massive loop 
contribution in calculations of interference effects. For the lightest Higgs mass considered 
here, (mu = 120 GeV) the prediction for the Higgs cross section including interference 
effects ((TH,i) is 10 — 15% lower than would be anticipated with the usual approach (<r#). 
Since a 0(10%) correction competes with the theoretical uncertainty of the NNLO cross 
section these pieces should be included in the prediction for a light Higgs cross section at 
the LHC. For mu > 400 GeV the interference becomes large and constructive, although 
one should take care in interpreting the results in this region since for these values of mu 
the Higgs width becomes very large (c.f. Table [|). For instance, including an s-dependent 
width for the Higgs boson would change the details of the results at large mu- 

A detailed examination of the interference can be performed by separating the Breit- 
Wigner form of the Higgs boson propagator into its real and imaginary parts. Thus the 
result for the interference can be written as, 

6(Ti = (s - m\f + m H T H *' r^s^ox} + ( ~ _ ^y + m H F H 3m \ 2 ^ m ^ A ^\ » 

(4.2) 
where -Affiggs = {s — m H + imu^u) Amggs- Although in the above equation the depen- 
dence of -Affiggs arid .Abox on the kinematic variables has been suppressed, both .Affiggs and 
Ahox contain thresholds and other kinematic structures. The results for the two terms in 
Eq. ( |4.2| ) are shown in Fig. |7| as a function of the invariant mass of the four final state 
leptons, m4£. Although mu is not an experimentally-measurable quantity, it is particu- 
larly important since it is equal to s, the partonic centre of mass energy, at this order. 
We show results for Higgs masses of 120 and 600 GeV. For mu = 120 GeV the width 
of the Higgs boson is very small (3.6 MeV) and hence the contribution of the imaginary 
part is negligible. The kinematic structure in the real part means that the s-dependence 
extracted in Eq. ( [4.2[ ) is not a faithful representation of the full dependence. Consequently, 
the dominant contribution to the interference comes from the real part of the Breit-Wigner 
as can be seen in Fig. |^. For mu = 600 GeV the imaginary part contributes about a 
third of the net interference. Although some cancellation between the region s < m H and 
s > m H is apparent, the net contribution is positive. For reference, in Fig. || we also show 
the contribution of the real part computed using two massless generations of quarks, which 
exhibits less kinematic structure. 

In Fig. H we present the distribution of mu, combining the Higgs signal (an or au,i) 
with the background (erg), for two different Higgs masses. We have chosen to compare the 
quantities d{au + &b) /dm^ and d(au,i + &b) /dm^ so that the distribution is positive 
throughout. In Refs. [13, 14] it was shown that, for the case of ZZ production, unitarity 
requires destructive interference between the Higgs and non-resonant contributions in the 
limit of large s. From our study of the interference as a function of mu, depicted in Fig. |6|, 
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Figure 7: Contributions to the interference coming from the real (red) and imaginary (blue) part 



of the Higgs Breit-Wigner propagator, as defined in Eq. (|4.2|), for two values of the Higgs mass. 
The real part of the interference with two generations of massless quarks is shown by the dashed 
magenta line. 



one might worry that the constructive interference at large mu violates the s requirement. 
However, Fig. |8| illustrates that this is not the case. Although there is a net constructive 
interference for run = 600 GeV, this is due to the large enhancement of the cross section 
in the region s < ran- In the limit of large s we observe that the interference is large and 
destructive, in complete agreement with the findings of refs. [13, 14]. For a much lighter 
Higgs boson, for instance ran = 120 GeV, there is no significant increase of the cross 
section in the region s < ran and the net destructive contribution simply results from the 
destructive interference in the tail. 
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Figure 8: The upper panels present results for the invariant mass of the four final state lep- 
tons using <jh + &b (blue dotted) and (7h,i + &b (red solid) for two different Higgs masses 
(left niH = 120, right mn = 600). In addition we plot the distribution with no interfer- 
ence terms but using the rescaled Higgs propagator, a is a + ob [10, 37] (black dashed). In the 
lower panels we present the ratio r = [d(an,i + o b) / dmu] / [d{(JH + aB)/dmi{] (red solid) and 
r = [d(a IS A + (J B )/dmu] / [d(a H + crB)/dmu] (dashed) . 



A method recently employed in ref. [10] to attempt to approximate the interference 
effects is based on the "improved s-channel approximation" (ISA) of ref. [37]. The prescrip- 
tion for unitarizing the amplitude is obtained by modifying the Higgs boson propagator, 



is 



im 



H 



m 



H s-m H +iT H {m H )^- 



(4.3) 



m H 



In order to test this approximation, in Fig. |§] we also show the m^ distribution calculated 
using this approach, denoted 07571- We observe that, as expected, the modified propagator 
decreases the cross section in the limit of large mu (s). We note that although the ISA 
approach has some of the features of the interference that we observe, there are significant 
differences in shape between the two approaches across the entire mu range. 

4.2 Interference with search cuts 

We now investigate the effect of the interference with more realistic search cuts mimicking 
those used by the CMS and ATLAS collaborations. To provide results for the CMS collab- 
oration we use the cuts that were employed in the search for the Higgs boson in the 2010 
data set [1]. Basic acceptance cuts are always applied to the missing transverse energy 
(rft) and lepton rapidities, 



tft > 20 GeV, \qt\ < 2.5 , 



(4.4) 



and then a number of further cuts are optimised for different values of the Higgs mass. In 
particular, cuts on the transverse momenta of the two leptons (pjJ iax , p^ m ), the invariant 
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m H [GeV] 


V h T [GeV] 


P S T [GeV] 


m cut [GeV] 


Ac^cut 


130 


25 


20 


45 


60° 


160 


30 


25 


50 


60° 


200 


40 


25 


90 


100° 



Table 6: Kinematic cuts used by CMS (taken from Ref. [1]) for a selection of potential Higgs 
masses. Cuts represent the lower limits on the lepton pr for the hardest and softest leptons, the 
upper limit on the invariant mass of the leptons mu and the maximum azimuthal angular separation 
between the leptons, A<pu (c.f. Eq. (4.5)). 



mass of the lepton pair (mu) and the azimuthal angle between the leptons (Acfiu) are all 
dependent on the value of mu that is assumed in the search. These cuts take the form, 



Pt > Pt > Pt > Pt 



m u < m cut , Acpu < A(f> c 



(4.5) 



where the cut thresholds for some benchmark values of ran are presented in Table ^. 

The cuts adopted by the ATLAS collaboration [3] are as follows. The set of basic 
acceptance cuts is, 



p lma* > 20 GeV, p l rf in > 15 GeV, tf t > 30 GeV, |%| < 2.5 



(4.6) 



and a further two cuts are applied that depend only on whether or not the putative Higgs 
boson is lighter than 170 GeV, 



m a < 50 (60) GeV, A<j) U < 1.3 (1.8) , 



(4.7) 



(numbers in parentheses indicate the values for mu > 170 GeV). In addition, the transverse 
mass Mt is constrained to be in the region 0.75 mn < Mt < mH, where Mt is defined as, 



M T 



(Eg + Ef iss ) 2 - (p*? + p^ iss ) 2 



(4.8) 



and E e T £ 



(p«)2 + 



m„ 



Our results for the cross sections with each set of cuts, with and without the interference 
terms, are depicted in Fig. || and tabulated for three values of ms in Table 0. For this study 
we have focused primarily on lighter Higgs bosons, where the H — > WW search channel is 
most relevant, and have considered three different choices of scale, (m#/2, m#, 2m#) when 
evaluating our predictions. At first sight the results appear qualitatively similar. As might 
be expected, the relative effect of the interference is insensitive to the choice of scale and 
grows with mn , as already observed for the case of no cuts in Fig. ^. Although the ratio 
of au,i to an is similar for both sets of cuts above 160 GeV, there are major differences 
in the overall impact of the interference terms in the low mass region, mu < 140 GeV. In 
this region one finds that, using the CMS cuts, the interference is destructive and leads to 
a reduction of around 10% in the Higgs signal cross section. On the other hand, the result 
of the interference terms after the application of the ATLAS cuts is an order of magnitude 
smaller, around 1%. 
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Figure 9: Cross sections for gg — >• H — >• W + (— > 
without interference effects using the CMS (left) and ATLAS (right) Higgs search cuts [1, 3] The 
red curves include the effects of the interference whilst the blue curves do not. Wc evaluate the 
cross section at three different scales m#/2 (dashed), run (solid) and 2m h (dotted). The lower 
panel shows the ratio between the cross sections with and without the interference, evaluated using 
each of the three scale choices. 



Experiment 


CMS 


ATLAS 


M H [GeV] 


130 


160 


200 


130 


160 


200 


&H 


1.58 


7.18 


3.02 


2.82 


9.35 


2.00 


<?H,i 


1.48 


7.16 


3.33 


2.79 


9.34 


2.24 


CTfT 


0.94 


1.00 


1.10 


0.99 


1.00 


1.12 



Table 7: Cross sections for gg — > H — > W + (— > v e e + )W~ (— > [i~v^) in femtobarns at i/s = 7 TeV 
with typical experimental Higgs search cuts applied (see Ref. [1, 3]), computed at leading order 
and either excluding (oh) or including {<JH,i) the effect of interference with the gluon-initiated 
background process. 



The large destructive interference effects under the CMS cuts are not surprising since 
similar results are found when no cuts are applied, c.f. Fig. |(J. We now investigate the 
ATLAS cuts in order to understand the origin of these destructive effects in the low mass 
region and why they are absent for this choice of cuts. The only major difference between 
the ATLAS and the CMS (or no) cuts is the application of the cut on the transverse mass 
of the system. We therefore present the Mt distribution for mu = 120 GeV, with no cuts 
applied, in Fig. 10. It is clear that the Mt distribution has a kinematic edge at Mt = ran 



since, if the Higgs were always produced on-shell, the transverse mass Mt would be zero 
at this order in the region Mt > van- However, small finite width effects ensure that this 
distribution has a small (but non-zero) tail. This argument holds for the resonant Higgs 
contributions to the amplitudes, but the non-resonant box diagrams have no such kinematic 
edge and instead can be relatively large in the Mt > run region. As a result the interference 
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Figure 10: The Mt distribution for a Higgs mass of 120 GeV with no cuts applied (upper panel). 
The lower (red) histogram is calculated including the effects of the interference term whilst the 
upper (blue) curve is obtained using the Higgs amplitude only. The ratio of the red and blue curves 
is shown in the lower panel. 



terms in this region have a large effect. In the region My > run we see that the prediction 
for what appears to be a cross section is actually negative. However ajj,i is not a true 
physical cross section, it is merely the change in the cross section induced by the presence 
of the Higgs boson in the calculation of the relevant amplitudes. To obtain a physical 
cross section, one would have to add the contribution from the background squared, ob- 
This contribution is large and positive and as a result the physical cross section remains 
positive. The destructive interference at large Mt is of course just a reflection of the 



expected destructive interference at large s already discussed in Section 4T and illustrated 
in Figs. ^ and ^. 

From this figure we can therefore conclude that the Mt cut employed by ATLAS 
naturally removes the region of phase space in which the interference is destructive and, 
as a result cth^/^h ~ 1 under the ATLAS cuts. It is also important, however, to quantify 
the effect of the interference on other kinematic distributions. One important example 
is the azimuthal angle between the leptons (A</>^) which is used in the search cuts of 
Eqs. (p~5,4.7). Since the Higgs signal peaks at low A<pi£ and the background peaks in the 
region close to A<pu = tt, this variable is a strong discriminant between signal-like and 
background-like events. Our predictions for this distribution, using ran = 120 GeV, are 



shown in Fig. |ll|. We observe that with no cuts applied there is a non-trivial shape change 
in A(p£ 



However, after the application of the cuts in Eqs. (|4.6, 4.7) and the Mt constraint, 
the shape change is negligible. This suggests that it is safe to use the ratio oh^I^h as an 
overall re-weighting constant under these cuts. 
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Figure 11: The A</>« distribution for a Higgs mass of 120 GeV, both with no cuts applied (left) 
and using the ATLAS cuts (right). The lower (red) histograms are calculated including the effects 
of the interference {<jh,%) while the upper (blue) curves are obtained using the Higgs amplitudes 
(oh) only. In the lower panel we plot the ratio of the two curves. 



To summarise, we would advocate the use of an upper limit on My in Higgs searches in 
the WW channel, for instance as currently used by the ATLAS collaboration. In addition 
to serving as a cut that reduces Standard Model backgrounds, it also reduces the destructive 
interference between the SM and Higgs amplitudes by around 10%, thereby increasing the 
expected Higgs signal. 



5. Higgs — y WW interference effects at the Tevatron 

The CDF and DO experiments at the Tevatron also use the WW channel to search for the 
Higgs boson in the mass range 120 — 200 GeV [2]. Since we are discussing gluon-gluon initial 
states the difference between the pp and pp colliders is immaterial and for this analysis the 
Tevatron and the LHC differ only in their centre of mass energies. The Higgs cross section 
at the Tevatron has been studied in great detail [4-10] , resulting in a theoretical uncertainty 
that is taken to be O(10%) [2]. Since the interference effects that we have studied so far 
at the LHC contribute at precisely this level, in this section we perform a more detailed 
study in order to quantify the effects in the Tevatron search region. 

We examine the effect of the interference both with no cuts and with typical search 
cuts used in the experimental analyses. Since the cuts used by CDF and DO do not 
differ significantly, we choose here to focus on the set of cuts employed in recent CDF 
analyses [2, 38]. Specifically, we require, 



p^ ax > 20 GeV, 



\rf max \ < 0.8, 

\rf min \ < 1.1, m u > 16 GeV 



(5-1) 
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Figure 12: Cross sections for gg — > H — > W + (— > v e e + )W~(-^ M _ ^m) m fcnitobarns at the 
Tevatron, with and without interference effects. The cross sections are computed by applying the 
Higgs search cuts (lower curves) or not (upper curves). The red curves include the effects of the 
interference whilst the blue curves do not. The lower panel shows the ratio cth^/ch for no cuts 
(solid), search cuts (dashed) and search cuts with an additional constraint Mt < mjj (dotted). 



In addition the missing transverse momentum is constrained using the ^ spec variable 
defined by [38], 



% 



spec 



^ 4 sin 



mm 



A<f>, 



TV 



(5.2) 



A(j) is the distance between the Ifc vector and the nearest lepton or jet. We require that 
$ 4 spec > 25 GeV. 



Our results both with and without this set of cuts are presented in Fig. 12. We observe 
that with no cuts the interference has the same structure as at the LHC. For mu < 2m\\r 
the effect is destructive and of the order 5% for mn < 130 GeV. Around run ~ 2my/ 
the effects are very small whilst for higher mu the effect is around 5% and constructive. 
We also observe that applying the search cuts increases the magnitude of the destructive 
interference in the low mass region. This effect, at the level of 10% for mu = 120 GeV, is 
as large as the theoretical uncertainty estimated from the NNLO cross section. The effect 
of the interference on current Tevatron analyses should therefore be taken into account. 

In the previous section we observed that a cut on the transverse mass Mt effectively 
eliminated the region of destructive interference. We have therefore recalculated the quan- 
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tities (Jh,i and an at the Tevatron, using the same CDF cuts described above but with an 
additional requirement that Mt < tuh- Under these cuts we indeed find that the impact 
of the interference is reduced considerably, from 0(10%) to 0(1%) just as at the LHC. 
The results for the interference with the additional Mt cut are illustrated by the dotted 



blue curve in Fig. 12 



6. Combining the interference with NNLO predictions 

In this paper we have computed the interference between Higgs and continuum contribu- 
tions at leading order, i.e. between the production of a Higgs boson via a top (or bottom) 
quark triangle and a gg —> WW box diagram. To calculate the interference at the next 
order one would need the NLO corrections to the Higgs process (which are known) and 
the NLO corrections to the gluon-gluon process (which are not). The NNLO result for the 
Higgs process has of course been known for some time [4-10], with large (factor of two) 
corrections to the LO cross section. 

A natural question arises as to how our results for the interference should be combined 
with the NNLO Higgs production cross section. One could simply modify the NNLO cross 
section by adding the LO interference terms, i.e. by defining, 

NNLO NNLO , t LO LO\ (R -, \ 

Since a^ NLO is much larger than ojj? one would expect that using Eq. ( |6.1| ) to estimate 
the effect of the interference would reduce the results in this paper by approximately a 
factor of two. 

However, treating the interference as an absolute correction to the NNLO cross sec- 
tion neglects the effect of all higher orders on the interference. Although the Higgs process 
receives large corrections at NLO and NNLO, the corresponding corrections to the gluon- 
gluon box diagram processes are unknown. Since we are not able to quantify these correc- 
tions a more conservative approach (particularly in the low-mass region where interference 
effects are large) would be to re-weight the NNLO cross section by the relative effect of the 
interference at LO, i.e. define, 

/rr LO \ 
NNLO NNLO H,i \ ,„ ^ 

It is this approach that we recommend in the region where the interference is destructive. 

7. Conclusions 

In this paper we studied the production of WW pairs through gluon fusion in the Standard 
Model. This process proceeds at loop- level through a closed fermion loop. In a previous 
paper [24] we presented analytic results for massless fermions, which correspond to the first 
two generations circulating in the loop. We extended these results to include contributions 
from the third generation keeping the full dependence on the top mass rrif. These formulae 
have been implemented into MCFMv6.1 which will soon be publicly available. The process 
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gg —)■ WW —)■ ul + l' v' keeping the full top (and bottom) mass dependence was previously 
calculated in ref. [12] using a semi-numerical approach. 

The amplitude for gg —> WW with non-zero mt can be written in the usual one-loop 
expansion in terms of scalar integrals. We found that the rational terms and all but two 
bubble coefficients were identical to the massless result. Further, the remaining two bubble 
coefficients were found to sum to the equivalent mt = coefficient. As a result we only 
needed to calculate one bubble coefficient. We used generalised unitarity techniques to 
calculate the box and triangle coefficients, which differ from the massless results. 

One of the most important phenomenological applications of these results is the effect 
on the Higgs/SM continuum interference. We illustrated that in general there are big 
deviations between the 2 and 3 generation calculations of this interference. It is crucial to 
quantify the impact of the interference as a function of the putative Higgs mass and under 
a variety of cuts, since in general these terms can contribute changes in rates comparable 
to the theoretical uncertainty associated with a NNLO prediction for the cross section. We 
found that in general the interference is largest far away from the 2mw threshold. For 
large m# (> 2mt) the interference is constructive and for low (m# < 140 GeV) masses it 
is destructive. In the region 120 < tuh < 130 GeV the destructive interference reduces the 
total cross section by around 10—15%. This large interference originates from non-resonant 
contributions above tuh- 

We illustrated that the interference is sensitive to the experimental cuts under consid- 
eration and that cuts used in the 2010 ATLAS analysis [3] dramatically reduce the large 
destructive interference for small ran- We showed that this is a result of the Mt cut that 
is employed, in particular the upper bound that effectively removes the non-resonant con- 
tribution to the interference. After applying the cut My < itih, the interference is reduced 
from O(10%) to 0(1%). Therefore we would advocate the use of the cut Mt < tuh for all 
experiments at hadron colliders since this enhances the signal by removing the region in 
which the destructive interference dominates. 

The total gg — >■ H cross section is known to NNLO [4-10] accuracy and the X-factor 
from going from LO to NNLO is large. Therefore, if the interference is treated as an absolute 
correction to the NNLO cross section its impact is reduced by an approximate factor of 
two. However, the interference is present at every order in perturbation theory and as such 
a more conservative approach would be to re- weight the NNLO cross section by the relative 
effect of the interference at LO. The size of the interference means that this is particularly 
important when setting Higgs mass limits using this channel for run ^ 140 GeV. 
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